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1 Introduction 

{Sj , First order PDEs of the form 

O 

^ = E^(«)^. i = l,...,n, (1.1) 

are called hydrodynamic or dispersionless systems in (l+l)-dimension. In this paper, 
^ we are interested in those PDEs among fll.lp which have multi-Hamiltonian structure, 



infinite hierarchy of symmetries and conservation laws. The important subclass of such 
systems are dispersionless limits of soliton equations. Differential Poisson structures for 
hydrodynamic systems were introduced for the first time by Dubrovin and Novikov [[J 
in the form 

7r y = ^(«)a,-X;r?(«)^ (i-2) 

k 

where g %3 is a contravariant flat metric and are related coefficients of the contravariant 
Levi-Civita connection. Then they were generalized by Mokhov and Ferapontov M to 
the nonlocal form 
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in the case of constant curvature c. The natural geometric setting of related bi-Hamiltonian 
structures (Poisson pencils) is the theory of Frobenious manifolds based on the geometry 
of pencils of contravariant metrics ||. 

The other methods of the construction of dispersionless systems are based on the 
application of the quasi-classical limit to the soliton theory. For example, the quasi- 
classical limit of dressing method is considered by Takasaki and Takebe |4[], while the 
quasi-classical limit of the scalar nonlocal <9-problem is presented by Konopelchenko and 
Alonso 0; see also the rich literature quoted in these papers. 

In the following, we develop an alternative approach to construction dispersionless 
systems and related Poisson pencils, based on an .R-matrix theory. As it is well known, 
the R- matrix formalism proved very fruitful in systematic construction of soliton systems 
(see for example |§]-[@] and the literature quoted therein). So, it seems reasonable to 
develop such a formalism for dispersionless systems. Recently, an important progress in 
that direction was made by Li [H who applied the .R-matrix theory to Poisson algebras 



10| . In this paper, we apply his results to a particular class of Poisson algebras. 

The paper is organized as follows. In section 2 we briefly present a number of 
basic facts and definitions concerning the formalism applied. In section 3 we apply the 
formalism of the classical .R-matrix to the Poisson algebras of formal Laurent series. 
Then in section 4 we illustrate our results with the known and new (l+l)-dimensional 
integrable dispersionless systems. 



2 Hamiltonian dynamics on Lie algebras: 
/^-Structures 

Let g be a Lie algebra, g* the dual algebra related to g by the duality map (•,•)—>• R, 
g* x g — >• R : (at, a) i— ► (a, a), a G g, a G g*, (2.1) 

and V(q*) := C°°(g*) be a space of C°°-functions on g*. Then, let 

ad : g x g — > g : (a, b) \— > ad a 6 = [a, b\, a, b G g, (2.2) 

be adjoint action of g on g, i.e. the Lie product, where [•, •] is a Lie bracket on g. There 
exists a natural Lie-Poisson bracket on T>(q*). Let F G X?(g*), then a map dF : g — ► g 
such that 

= (L',dF(L)), L,L'e & *, (2.3) 

t=o 

is a gradient of F. Let L G g*, functions H,F belong to the space of functions on 
g* : T>(q*), and their gradients dH, dF G g, then the Lie-Poisson bracket reads 

{H,F}(L) := (L, [dF, dH}). (2.4) 

We confine our further considerations to such algebras g for which its dual g* can 
be identified with g. So, we assume the existence of a product (•, -) on g which is 
symmetric, non-degenerate and ad- invariant: 



-F(L + W) 



(ad a 6, c) + (6, ad a c) = 0, a,b,ce$. 



(2.5) 
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Then, we can identify g* with g, (g* = g) by setting 

(a, 6) = (a,6) fls a,beg,aeg*, (2.6) 

where a G g* is identified with a G g. Now, we can write the Lie-Poisson bracket as 

{H,F}(L) = (L, [dF,dH]) = (L, [dF,dH]) e 

= (dF, [dH, L]) a = (dF, -ad L dH) s = (dF, 6(L)dH) s , (2.7) 

where 9 is a Poisson tensor 9 : g — > g*. Hence, the Hamiltonian dynamical system on 
g* can be defined by the equation 

L t = 9(L)dH = -ad L dH = [dH, L). (2.8) 



Now, we can identify the dynamic equation (|2.8 ) and the Lax equation with a natural 
Hamiltonian structure 

L t = [A, L] = 9(L)dH = [dH, L]. (2.9) 

This abstract approach to integrable systems profits from a deeper understanding of 
the nature of integrability as well as equips us with a very general and efficient tool for 
the construction of multi-Hamiltonian systems from scratch. 

Definition 2.1 An R-structure is a Lie algebra g equipped with a linear map 
R '■ Q ~^ (called the classical R-matrix) such that the bracket 

[a,b] R := [Ra,b] + [a,Rb], a, b G g, (2.10) 

is a second Lie product on g. 

Definition 2.2 Let A be a commutative, associative algebra with unit 1. If there is a 
Lie bracket on A such that for each element a G A, the operator ad a : b i— > [a, b] is a 
derivation of the multiplication, then (A, [.,.]) is called a Poisson algebra. 

Thus the Poisson algebras are Lie algebras with an additional associative algebra struc- 
ture (with commutative multiplication and unit 1) related by the derivation property 
to the Lie bracket. 

Theorem 2.3 j^j Let A be a Poisson algebra with Lie bracket [-, ■] and non- degenerate 
ad-invariant pairing (-,-)a with respect to which the operation of multiplication is sym- 
metric, i.e. (ab,c)A = (a, be) a, Va, b, c G A. Assume R G End(A) is a classical R- 
matrix, then for each integer n ^ — 1, the formula 

{H, F} n = (L, [R(L n+1 dF), dH] + [dF, R{L n+1 dH)}) A , (2.11) 

where H, F are smooth functions on A, defines a Poisson structure on A. Moreover, all 
{-,-} n o,re compatible. 
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The related Poisson bivectors 7r n are given by the following Poisson maps 

ir n :dH^-ad L R(L n+1 dH)-L n+1 R*(ad L dH), n ^ -1, (2.12) 

where the adjoint of R is defined by the relation 

(a,Rb) A = (R*a,b) A . (2.13) 

Notice that the bracket ( |2.11| ) with n = — 1 is just a Lie- Poisson bracket with respect 
to Lie bracket ( |2.10| ) 

{H,F}_ 1 = (L,[dF,dH] R ) A . (2.14) 

We will look for a natural set of functions in involution w.r.t. the Poisson brackets 
(|2.11|) . A smooth function F on A is ad-invariant if dF e ker adL, i.e 

[dF,L]=0, LeA, (2.15) 

which are Casimir functionals of the natural Lie- Poisson bracket ( |2.4j) . 
Hence, the following Lemma is valid 

Lemma 2.4 j^j Smooth functions on A which are ad-invariant commute in {•, -} n . The 
Hamiltonian system generated by a smooth ad-invariant function C(L) and the Poisson 
structure {•, ■}„ is given by the Lax equation 

L t = [R{L n+1 dC),L], LeA. (2.16) 

Let us assume that an appropriate product on Poisson algebra A is given by the 
trace form tr : A — > R 

(a,b) A = tr(ab), a,b G A. (2.17) 

As we have assumed a nondegenerate trace form tr on A, we will consider the most 
natural Casimir functionals given by the trace of powers of L, i.e. 

C q (L) = -L-tr(L« +1 ). (2.18) 



The related gradients by ( |2.3| ) are of the form 

dC q (L) = L q . (2.19) 

Then taking these C q (L) as Hamiltonian functions, one finds a hierarchy of evolution 
equations which are multi-Hamiltonian dynamical systems 

L tq = [R(dC q ), L] = 7r_i(d<7,) = 7ro(dC 5 _i) = ... = ^(dC^) = ... . (2.20) 

For any i?-matrix each two evolution equations in the hierarchy ( |2.20 ) commute due 
to the involutivity of the Casimir functions C q . Each equation admits all the Casimir 
functions as a set of conserved quantities in involution. In this sense we will regard 
Q2.20D as a hierarchy of integrable evolution equations. 
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To construct the simplest /^-structure let us assume that the Poisson algebra A can 
be split into a direct sum of Lie subalgebras A + and v4_, i.e. 

A = A + ®A_, [4 ± ,4]c4 (2.21) 

Denoting the projections onto these subalgebras by P±, we define the /2-matrix as 

R=\{P + -P-) (2.22) 

which is well defined. 

Following the above scheme, we are able to construct in a systematic way integrable 
multi-Hamiltonian dispersionless systems, with infinite hierarchy of involutive constants 
of motion and infinite hierarchy of related commuting symmetries, ones we fix a Poisson 
algebra. 



3 Poisson algebras of formal Laurent series 



Let A be an algebra of Laurent series with respect top [11 



I i& J 

where the coefficients Ui(x) are smooth functions. It is obviously commutative and 
associative algebra under multiplication. The Lie-bracket can be introduced in infinitely 
many ways as 

[Lu L a ]=p{——- —— j := {L U L 2 } r , r e Z, (3.2) 

as ad^ = P r (^-§^ — f^rj^) is a derivation of the multiplication, so A r := (A, {•, -} r ) are 
Poisson algebras. 

Lemma 3.1 An appropriate symmetric product on A r is given by a trace form (a, o)a '■ = 
tr(afe): 

trL = / res r Ldx, res r L = u r _i(x), (3-3) 



which is ad-invariant. In expression (|3.3j ) the integration denotes the equivalence class 
of differential expressions modulo total derivatives. 

Proof. We assume that Q = S 1 if u is periodic or Q = R if u belong to the Schwartz 
space. The Symmetry is obvious as LiL 2 = L 2 L X . Let L 1; L 2 £ A : L\ = aip\ L 2 = 
Y^j b j p j , then 



res r [L u L 2 ) = res r lp r ^2{ia i {b j ) x -j(a i ) x b j )p l+: > 1 J = ^{(ajb. 



(3.4) 
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So, tr[L 1; L 2 ] = and hence 



tr{[A,B]C)+tr{B[A,C}) 

= tr([A, BC] - B[A, C}) + tr(B[A, C}) = tr[A, BC] =0. □ 
For a given functional F(L) = f Q f(u)dx, we define its gradient as 

i 

where Sf/Sui is a variational derivative. 

We construct the simplest i?-matrix, through a decomposition of A into a direct sum 
of Lie subalgebras. For a fixed r let 



A^- r+ k = P^-r+kA = < L = ^2 u i( x )p l \ 

A<- T+ k = P<- r +kA = I L = ^ U i( X )p r \ 

I i<-r+k ) 



(3.6) 



where P are appropriate projections. 

Proposition 3.2 A^_ r+ k, A < _ r+ k are Lie subalgebras in the following cases: 

1. k = 0, r = 0, 

2. k = 1,2, r G Z. 

The proof is through a simple inspection. Then, the i?-matrix is given by the projections 

R — ~^(P^-r+k — P<-r+k) = P^-r+k ~ 2 = 2 ~ ^ >< ~ r+k ' (3-7) 

To find R* one has to find P^_ r+k and P<_ r+fc given by the orthogonality relations 

(P^-r+k: P<-r+k) = (-P<_ r+fc , P^-r+k) = 0. (3.8) 

So, we have 

R^-r+k — R<2r-k, R<-r+k = R^2r-k, (3.9) 

and then 

R* = 2~(R>-r+k - R<-r+k) = 2 ~~ R>2r-k = P<2r-k ~ ^ (3.10) 

Hence, the hierarchy of evolution equations ( |2.20| ) for Casimir functionals C (L) with 
i?-matrix given by (|3.7|) has the form of two equivalent representations 

L tq = {(L%„ r+fc , L} r = -{(L <? ) < _ r+fc , L} r , L e A, (3.11) 
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which are Lax hierarchies. 

We have to explain what type of Lax operators can be used in (|3.11| ) to obtain a 
consistent operator evolution equivalent with some nonlinear integrable equation. Here, 
we are interested in extracting closed systems for a finite number of fields. The case of 
infinite number of fields was considered recently in [IT]. Hence, we start with looking 
for Lax operators L in the general form 



unP N + Wat-iP 1 + ... + U- 



^p' m+1 + u. 



(3.12) 



of N-th order, parametrized by finite number of fields Uj. To obtain a consistent Lax 
equation, the Lax operator (|3.12|) has to form a proper submanifold of the full Poisson 
algebra under consideration, i.e. the left and right-hand sides of expression ( |3.11| ) have 
to lie inside this submanifold. 

Observing (|3.11| ) with some (L 9 )<_ r+ fc = a_ r -ffc_ip~ r+fe_1 + a^ r+k -2P~' +k ~ 2 + ■■■ one 
immediately obtains the highest order of the right-hand side of Lax equation as 

L t = (u N ) t p N + (u N ^) t p N ' 1 + ... 



-{(L 9 )<_ r+fc , u N p N + lower}, 



-((-r + k - l)a- r+k -i(u N ) x - N(a- r+k -i) x u N )p 



N+k-2 



+ lower, 



(3.13) 



where lower represents lower orders. Observing ( |3.11| ) with some (L q 



+ 



a-r+k+iP r+fc+1 + a_ 
side of Lax equation ( |3.11| ) as 

u 



+kP r one immediately obtains the lowest order of the right-hand 



+ (u. m+1 ) t p- m+1 + (u„ m ) t p- m 



{(L q 



<^-r+k 



, higher + u_ m p m } r 



= higher + ((-r + k)a_ r+k (u_ m ) x - {-m)(a_ r+k ) x u^ m )p m+k l , (3.14) 

where higher represents higher orders. Simple consideration of ( |3.13|) and (|3.14|) with 
the condition iV ^ — m leads to the admissible Lax polynomials with a finite number 
of field coordinates, which form proper submanifolds of Poisson subalgebras. They are 
given in the form 

(3.15) 
(3.16) 
p~ m , (3.17) 

are arbitrary time independent 



k 
k 
k 



0, r = 

1, r e Z 

2, reZ 



L 
L 



C N p N + C N ^p N 1 + U N ^ 2 P N 2 + ■■■+ U x p + M , 
cnP N + u N -ip N ~ l + ... + U\- m p l ~ m + u- m p~ m , 



unP N + un-iP N 1 + ... + Ml 



where the Ui are dynamical fields and cjv,cjv-i,C- 
functions of x. 

Once we know the restricted Lax operators L G A, we can now investigate the form 
of gradients of Casimir functionals given by powers of L, as well as we can investigate 
some further simplest admissible reductions of Lax operators. 

In general, the fractional powers of L will lead to interesting results. Let L be given 
by ( p. 12 ), then we consider polynomials of the form 



-i 



a x p + a + a_ip 
a-ip" 1 + a- 2 p~ 2 + a-3P~ 3 + 



... + bip + b 
... + hp 3 + b 2 p 



b- lp -\ 
2 ' hp, 



for N eZ^ 
for JVeZ- 
for m G 
for m G Z_ 



(3.18) 
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where the coefficients a, and b { are obtained by requiring (L~^) N = L and (L») ra = L, 
successively via the recurrent procedure. Therefore, one finds the formal expansion of 
Ln and I» and one can calculate the fractional powers of L for integer q: Ln and L™. 
Notice that they are in the form of infinite series, except the case of integer powers, 
obviously. In fact, we need only their finite parts (L^)^_ r+ fc or (Lm)<_ r+ fc. Hence, for 
a given L E A in principle we can construct two different hierarchies of Lax equations 



The case of k = 0. Let us consider Lax operators of the form ( |3.15| ). One can see 
that Ln has the form 

= a q p q + a q ^p q ' x + a q - 2 p g ' 2 + a q ^p q ~ 3 + lower, q E Z+, (3.19) 

where oti, are arbitrary x independent functions. The second form L™, since m — 0, 
gives only the integer powers of L, such that (£ 9 )^o — L q , leading to trivial dynamics 
L t = {L q , L} = 0. Hence, for k = there is only one Lax hierarchy for gradients of 
Casimir functionals ( |3.19| ). There are no further reductions. 



The case of k — 1. Let us consider Lax operators of the form ( |3.16| ). One can see 
that Ln and L™ have the forms 

L% = a q p q + a q ^p q ' 1 + a q . 2 p q ~ 2 + a q _ 3 p q ~ 3 + +lower ?6Z + , (3.20) 

Lm = higher + a^ q p i ~ q + a 2 - q p 2 ~ q + a^p 1 ^ + u™ m p~ q q G Z+, (3.21) 

where cti is an arbitrary ^-independent function. We remark that there is always a 
further admissible reduction of equations (|3.11|) given by u_ m = 0, since such reduced 
Lax polynomial would still be of the form (|3.16|). After such reduction we have to look for 



the form of gradients of Casimir functionals. By putting w_ m = in ( |3.20| ), it preserves 



the order of highest terms and the form. For ( |3.21|) the lowest order disappears, and as 
all other terms depend linearly on the powers of w_ m , such an L~ will reduce to zero, 
except the one case for q = m. We can continue the reductions by putting u\- m = and 
so on. Therefore, the reductions are proper in general only for the gradients of Casimir 
functionals in the form ( |3.20| ). 



The case of k = 2. Let us consider Lax operators of the form ( |3.17| ). One can see 
that Ln and I™ take the form 

= ufp q + a q „ x p q - x + a q . 2 p q ~ 2 + a q _ 3 p q - 3 + lower q E Z+, (3.22) 

Lm = higher + a 3 ^ q p 3 ~ q + a 2 ^ q p 2 ~ q + ax^ q p l ^ q + a- m p~ q q E Z+, (3.23) 

where is an arbitrary x-independent function. We remark that there is always a 
further admissible reduction of equations (|3.11|) given by = 0, since such a reduced 



Lax polynomials would still be of the form ( p,17|) . The next reduction is ujv-i = and 



so on. By analogous considerations as for k — 1, these reductions are proper in general 
only for the gradients of Casimir functionals in the form (|3.23| ). 

The different schemes are interrelated as it is explained in the following theorem. 
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Theorem 3.3 Under the transformation 

x' = x, p' = p~\ t! = t (3.24) 

the Lax hierarchy defined by k = 1, r and L transforms into the Lax hierarchy defined 
by k = 2, r' = 2 — r and V , i.e. 

k = l,r,L k = 2,r' = 2 - r,L'. (3.25) 



Proof. It is readily seen that the Lax operators for k = 1 and r of the forms ( |3.16| ) 
transform into the well restricted Lax operators for k = 2 and r' = 2 — r of the forms 
(|3.17|) . Let's observe that 

{a, B} r = - ^^) = -y- r+2 (— — - — — ) = ~{A', B'}'„ 

r dp dx dx dp dp' dx' dx' dp' 

and 

(dcy >s = (dc f u„ s . 

Hence, we have 

U = {(dC)^ r+1 ,L} r = -{(dC)'^_ r+1 ,L'}' r , 

= -{(dC% r ^,L'y r , = -{(dC'U^+^L'}',, = L' t ,. □ 



Therefore, some dispersionless systems can be reconstructed from different Poisson 
algebras. Moreover, we remark that the gradients of Casimir functionals for k — 1 of 
the form (|3.20|) , (|3.21| ) by p' = p^ 1 transform into ( |3.23|) , ( ^.22|) for k = 2, respectively, 
at a slant. 



Two equivalent representations of Poisson bivectors (p.l2j) with the i?-matrix given 
by ( p.7|) are defined through the following Poisson maps 

7r n dH = {(L n+1 dH)^ r+k , L} r + L n+1 ({L, dH} r )> 2r _ k 

= -{(L" +1 ^)<_ r+fc , L} r - L n+1 ({L, dH} r ) <2r _ k , n > -1. (3.26) 

It turns out that the first representation yields a direct access to the lowest polynomial 
order of ir n dH, whereas the second representation yields information about the highest 
orders present. There are two options. The best situation is when a given Lax operator 
forms a proper submanifold of the full Poisson algebra, i.e. the image of the Poisson 
operator ir n lies in the space tangent to this submanifold for each element. If this is not 
the case, the Dirac reduction can be invoked for restriction of a given Poisson tensor to 
a suitable submanifold. 

The case of k = 0. Let us consider the simplest admissible Lax polynomial ( |3.15| ) of 
the form 

L = p N + u N ^ 2 p N ~ 2 + ...+u lP + u . (3.27) 
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This is the well-known dispersionless Gelfand-Dickey case. Then, the gradient of the 
functional H(L) is given in the form 

SH SH i SH 2 SH , N 

TT = 1~P + 1~P + - + 1 P ■ ( 3 - 28 ) 

dL OUq OU\ OUn-2 



Observing ( |3.26| ) for n — — 1 one immediately obtains the highest and lowest order 
of Ti^idH as 

7T_i f = a N - 2 p N ~ 2 + ajv- 3 p JV " 3 + ... + Oi lP + a . (3.29) 

Hence ir_idH is tangent to the submanifold formed by the Lax operator of the form 
( |3.27| ), and the linear Poisson structure, since (ff-) >0 = 0, is given by 



SL J VI ' SL J 



o/ >o 



For n = 0, L does not define a proper Poisson submanifold, as 



SH 

Jl 



n o I tt" ) = "w-iP^ 1 + ®n-2P N 2 + ••• + otip + a , 



and a Dirac reduction is required. Let 

L = p N + up N ' 1 + u N „ 2 p N - 2 + ...+u = L + up"' 1 (3.31) 

be a an extended Lax polynomial and we shall consider the 7r Hamiltonian flow for 
L together with the constraint u — 0. However, imposition of such a constraint leads 
to the modification of the tc Poisson structure due to the Dirac reduction. We briefly 
remind the calculation procedure on the example considered. The Hamiltonian flow 
for u, given by the coefficient of p N_1 in the Hamiltonian equation for L, under the 
constraint u = 0, gives the relation 



^L=o= (^o{S,l|J =0, (3.32) 

where 

SH 



Then solving ( |3.32| ) with respect to j^- one gets 



sh i ! r 



(3.34) 
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It means that the function ^ can be expressed in terms of J^. This implies 



l S JL + s JL p -« 

' SL Su 



SH\ f8H\ \ S8H\ A ([ 5H 




0/ >o 



0/ >0 



>0 



^}J,4{^{^}.}„ 



(3.35) 

i.e. the second Poisson map of dispersionless Gelfand-Dickey systems. Poisson structures 
7r_! and 7TQ ed were constructed for the first time in |12| as the dispersionless limit of the 



Poisson structures of the Gelfand-Dickey soliton systems. Notice that ^ d is purely 
differential due to the property (|3.4j ). 

Observing Q3.26 ) for n ^ 1 one obtains the highest and lowest order of ir n dH as 



7r n (jpj = a^N-ipW- 1 + a (n+1)N _ 2 p^ N - 2 + ... + a lP + a . (3.36) 



Hence, the polynomials of the form ( |3.27|) do not form a proper Poisson submanifold. In 



fact there is not obvious proper Poisson submanifold for 7r n with n ^ 1, apart from the 
trivial case of the first order polynomials with n = 1. Nevertheless, the Dirac reduction 
can be invoked to restrict the bivectors n n on the polynomials to the submanifold of the 
form ( gggp . 



The case of k — 1. This case contains new results. Let us consider the simplest 
admissible Lax polynomial fl3.16| ) of the form 



L = p N + u N ^p N ~ l + ... + Ul . mP l - m + u. m p~ m . (3.37) 
Then gradient of the functional H(L) is given in the form 

5 4- = ^f^ 1 + -^ V ™ + + ~7~~ — p r ~ N ■ (3.38) 

OL 0U- m 0U- m+ i OUtf-l 



Observing (|3.26|) for n = — 1 one obtains the highest and lowest order of 7r_! (j^) as 



71". 



_ x ( 5 Jf) = ((...y^ 1 + lower) + ((...)p 2r - 2 + lower) 



\5L ) 



(higher + (...)p~ m ) + (higher + (...)p 2r_1 ) 



where lower (higher) represents lower (higher) orders. Hence, the Lax operators of 
the type ( |3.37| ) form a proper submanifold for N ^ 2r — 1 ^ — m, as then 7r_i (j£) is 
tangent to this submanifold. So the linear Poisson map is 

^{^) = {[^\ ,L> +[{L,^} ) . (3.39) 




r/ >2r-l 
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Otherwise a Dirac reduction is required. 

For the second Poisson map with n = 0, L does not define a proper Poisson sub- 
manifold and two distinct cases have to be considered. 

2r - 1 > 1 : 

tto (jp) = (...) P ^ + ^ + ... + (..Op"" 1 + ... + 

hence L is not properly defined and a Dirac reduction is required for additional higher 
order terms. The simplest case is r = 1 with one-field reduction. Let 

L = u P N + u N ^ lP N - 1 + u N ^ lP N - 2 + ... + «i-„y- m + u. m p' m . 

The Dirac reduction with the constraint u — 1 leads to the second Poisson map in the 
form 



7T, 



^red 




5# 



which is purely differential. 
2r - 1 < : 

tto 0|) = (..O^ 1 + - + (-) P - m + ••■ + (..>~ m+(2r ~ 1) , 

hence L is not properly defined and a Dirac reduction is required for additional lower 
order terms. The simplest case is r = with one-field reduction. Let 

L = p N + Un-xp"- 1 + ... + u^ mP l - m + u^ m p~ m + up~ m -\ 

The Dirac reduction with the constraint u = lead to the second Poisson map in the 
form 



which is purely differential. This special case was considered recently in [JT3[]. 

Observing Q3.26D for n ^ 1 one obtains the highest and lowest order of 7r n (j^) as 




^ ( 5 Jf) = ((...)p N -i + Wer) + ((...)p^ N+2r - 2 + lower) 



{higher + {...)p~ m ) + (higher + {...)p~ {n+l)m+2r - 1 ) 



where lower {higher) represents lower (higher) orders. Hence, the Lax operators of the 
type (|3.37 ) do not form a proper Poisson submanifold for the 7r n with n ^ 1, apart from 



the trivial case of N = —m = — ^. Hence, one has to apply Dirac reduction to restrict 



the bivectors n n on the polynomials to the submanifold of the form (|3.26| ) 
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The case of k = 2. This has not been considered yet. Let us consider a Lax polyno- 
mial (|3.17|) of the form 

L = u N p N + un^- 1 + ... + u^ mV x - m + p~ m . (3.42) 

Then, gradient of functional H(L) is given in the form 

6 4 = 4^V™ + ••• + -^f- N + ^f-"- 1 - (3.43) 

dL OUi- m OUN-1 OUn 



Then by analogous consideration as for k = 1 or by Theorem |3.3| , for the first Poisson 
structure with n = —1, L defines a proper Poisson submanifold for N > 2r — 3 > —m, 
so the first Poisson map in this case is 

-i ( 6 4) = UiT) L) +(U-7T-N • '3.11) 



\5L J 




r/ >2r-2 



Otherwise a Dirac reduction is required. 

For the second Poisson map with n = 0, L does not define a proper Poisson sub- 
manifold and again two distinct cases have to be considered. 



2r - 3 > : 

tto 0§) = (...)p n+ ^ + ... + (...) P N + ... + (...y- m , 

hence L is not properly defined and a Dirac reduction is required for additional higher 
order terms. The simplest case is r = 2 with one-field reduction. Let 

I = U p N+1 + U NP N + U^p"- 1 + ... + U^ mV X ~ m + p~ m , 

then the Dirac reduction with the constraint u = leads to the second Poisson map in 
the form 



„red 



SH\ 
~5LJ 



2 

which is purely differential 



, (3.45) 

2J 2 



2r - 3 < : 

tto (jfj = (...)P N + ... + (-)p- m+1 + ... + (...K m+1+(2 '- 3 \ 

hence L is not properly defined and a Dirac reduction is required for additional lower 
order terms. The simplest case is r = 1 with one-field reduction. Let 

L = U N p N + UjV-li^" 1 + . - + U 2 -mP 2 ~ m + U^rnp 1 '" 1 + Up' m , 
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then the Dirac reduction with the constraint u = 1 leads to the second Poisson map in 
the form 



'red 
7Tn 



6H_ 
~6L 



which is again purely differential. 

Now we present one example of three-field Dirac reduction. Let us consider the case 
with r = 0, then 

I = u N p N + ... + u^ mP l - m + up- m + vp- m - 1 + wp- m ~ 2 . 

The Dirac reduction with constraints u — 1, v = w = gives the following reduced 
Poisson map 

'6H\ 

(L S JL ) ^ , A + L ({ L, <* } J + { L, A P + B + Cp-X , (3.47) 



„red 
^0 



where 



5/n 



m VI J t 
5= V (7^11 



' 61 KJ -1 

If 1 TO — 1 2 \ 1 / 1 TO — 1 2 \ 

to \ 2 m J m \ 2 m J 

generally nonlocal. 

Then by analogous consideration as for k = 1 or by Theorem |3.3| , we see that, Lax 
operators of the form ( |3.42j ) do not form a proper Poisson submanifold for the Tc n with 
n ^ 1, apart from the trivial case of N = — m = ^z^n Hence, one has to apply the 
Dirac reduction to restrict the bivectors 7r n on the polynomials to the submanifold of 
the form (|3~26|) . 

Hence we know the Poisson structure for (l+l)-dispersionless systems constructed 
from Poisson algebras, and since we are interested in multi-Hamiltonian systems 

L tq = {(L«);>_ r+fc , L} r = n_ x dH x = n dH = 7r_idff_i = ... , (3.48) 

we shall now consider the problem of their construction. The conserved quantities Hi 
from (12.181) are defined as follows 



HAL) = -^—tr(L q+i ) = — *— I res r (L q+l ) dx. (3.49) 
q + i q + iJn 
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4 A list of some (l+l)-dimensional dispersionless 
systems 

In this section we will display a list of the simplest nonlinear dispersionless integrable 
systems. Calculating the powers L» we consider the Lax hierarchy 

L tn = [ {L%)^_ r+k ,l} r , n = 1, 2, 3, ... . (4.1) 

The second hierarchy with powers L™, can be obtained by the transformation from 
Theorem \5.'d[ which we leave for the interested reader. In general for simplicity we 
present only the bi-Hamiltonian structure. For k = and k = 1 the choice n = 1 — r 
will always lead to the dynamics (ui) tl _ r = (1 — r){ui) x for the fields Ui in L, so that 
we may identify t\- r = jz^x in this cases. For k = and integer values of n/N the 
equations become trivial, because then (L"^)^o — L. For each choice of k = 0, 1 or 2 
and N we will exhibit the first nontrivial of the nonlinear Lax equations (|4.1|) associated 
with a chosen operator L. 



The case of k = 0. 

Example 4.1 Dispersionless Korteweg-de Vries: k = 0, r = 0, iV = 2. 

This is a standard case of the dispersionless Korteweg-de Vries (dKdV) hierarchy. 
The Lax operator for the dKdV has the form 

L = p 2 + u. (4.2) 

We derive the dKdV equation 



J t 3 



| (L§) >q , L j ^ Mfe = ^ = TT-idHi = n r ed dH = ir^dH^, (4.3) 



where we get the Poisson tensors from ( |3.30| ) and (|3.35| ) 
7r_ 1 = 2a BJ vr red = d x u + ud x , 



7T 



" - n r ed (n^y 1 n r ed = d x u 2 + u 2 d x - Ui x d~ x u x , 



(4.4) 



2 

and the respective Hamiltonians 

H\ — - j v? dx, H Q = - j u 2 dx, H-x = f u dx. (4.5) 

Example 4.2 Dispersionless Bousinesq: k = 0, r = 0, iV = 3. 
The Lax operator is given by 

L = p 3 + up + v . (4.6) 

We derive 



£0 in V V 



Li) ,L\ =( ) = *-idH x = ^ d dH,. (4.7) 



3 



16 



Maciej Blaszak and Blazej M. Szablikowski 



Eliminating the field v from this equation we can derive the (l+l)-dimensional disper- 
sionless Boussinesq equation 
2 

u tt = -^(u 2 ) xx . (4.8) 
The respective Poisson tensors are 

o ( d x \ red ( d x u + ud x 2d x v + vd x \ , . 

n - 1 = 3 {d x J' *° -\d m v + 2vd g -\ud x u J' ^ 

and the Hamiltonians are given in the following form 

H x = - f(v 2 - -u 3 ) dx, H = I v dx. (4.10) 



3Jn 9 Jn 

Example 4.3 The three field case: k = 0, r = 0, N = 4. 
The Lax operator is 

L = p 4 + up 2 + vp + w, (4-11) 




K_ x dH x = 7T r dH , (4.12) 



where 



Ad x 
7T_i = | Ad x 

4d x d x u + ud x 

d x u + ud x 2d x v + vd x 3d x w + wd x 

iTn ed = | d x v + 2vd x —ud x u + 2d x w + 2wd x —\ud x v 

d x w + 3wd x —\vd x u —\vd x v + ud x w + wd x u 



The case of k = 1. 



Example 4.4 Three field hierarchy: k = 1, r E Z \ {2}. 

The Lax operator has the form ( |3.16|) with N = 2 — r, m = r + l 







• 






<2-r 





(4.13) 



#i = - / (-~it 2 f + fw) rfx, H = [ v dx. (4.14) 

2 4 y n 



L = p 2 " r + Mp 1 " r + vp~ r + u7T r " 1 , (4.15) 
Then we find 

(2 - r)^ 

r«xf + (1 - + (2 - r)w x | , (4.16) 
(1 + r)u x w + (1 — r)uw x 
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This Lax operator forms a proper submanifold as regards 7r_i only for r = 0, 1. Other- 
wise a Dirac reduction is required. Then for r = 

2v x 

uv x + 2w x \ = vr.xcffii = ir r ed dH , (4.17) 
u x w + uw 3 



V 






'.-( 


w 





where 



2a 

7T_! = I 2<9 X wo» 






49 x u 2<9. t .t> 

^red _ j 4:iid x 2ud x u + d x v + f a x U(9 x f + 2d x w + wd x 
vd x u + d x w + 2w(9 x + wd x u 



(4.18) 



H\ = vw dx, H = w dx. (4.19) 
Jn Jn 

For r = 1 we have 



V 






'.-( 


w 









where 




TT-idfTa = n r ed dH , (4.20) 



c^f 

7r_i = | fc^ c^w + wd x 

2^^ 

d x v + f (9 X + 2d x w + itfd z 2ud x w 

7To erf = ( f 9 2 m + 9jj.ii; + 2^(9^, 2^9^.^ + w^iu + wd x u Avd x w 

2wd x u Awd x v Qwd x w 



(4.21) 



Hi = ^ I {u 2 + 2v) dx, H = I u dx. (4.22) 

Example 4.5 Dispersionless Toda: k — 1, r G Z \ {2} 

The first admissible reduction w; = of ( |4.15| ) leads to the two field Lax operator 

L = p 2 ~ r + up 1 -" + vp~ r . (4.23) 

This Lax operator forms a proper submanifold as regards ir_i only for r = 1, in other 
cases a Dirac reduction is required. For r = 1 by reduction of ( |4.20| ) we get the disper- 
sionless Toda equation 

( " ) = ( Z ) = w -" m = <d " H " (4 - 24) 
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where 



d x v \ d _ f d x v + vd x ud x v 

\ ri) x ) ' V vd x u 2vd x v 

1 - .2 



Hi — - J (u + 2t>) dx, H = J u dx. 



n 



For r = we have 



u \ ( 2v x 
v J t2 \uv x 



(4.25) 



(4.26) 



but we lose the bi-Hamiltonian structure since there are no Dirac reductions with the 
constraint w = of ( |4.18|) . 



The next admissible reduction w = v = of ( 4.16|) leads to the noninteresting trivial 
equation L t2 r = since (L)^_ r+1 = L. 

Example 4.6 Three field hierarchy: fc = l,rGZ\{l}. 

The Lax operator has the form ( |3.16| ) with iV = 1 — r, m = r + 2 



L = + up~ r + vv~ r ~ x + wp~ r ~ 2 , (4.27) 
Then we find 



L t2 _ r = <! [ ) . /, 



2-r 



1 — r 

t'2-r 



UU X + (1 - r)^ 

1 + r)u x v + (1 - r)uv x + (1 - r> x J. (4.28) 
(2 + r)u x w + (1 — r)uw x 



This Lax operator forms a proper submanifold as regards 7r_! only for r = 0, in other 
cases a Dirac reduction is required. Then for r = we have 







• 




w 





MMi + V x 



U X V + UV X + W X I = 7T_ 1 rfi7 1 = TT^dHo, (4.29) 

2M3.it; + ttu; x 



where 






rr_i = i <),. 

c^-u; + w9 x 

fa, d x u \d x v 
n o ed = I u ^x 9 x v + vc^ 2S3.1t; + wd x 

\vd x d x w + 2wd x —\vd x u + ud x w + wd x u 



(4.30) 



H 1 = (u 2 v + v 2 + 2uw) dx, H = (uv + w) dx. (4.31) 
Jo. Jn 
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Example 4.7 Benney system: k = l,r G Z \ {1} 

The first admissible reduction w = of (|4.27|) leads to the two field Lax operator 

L = p 1 ^ + up~ r + vp' r -\ (4.32) 

This Lax operator forms a proper submanifold as regards n_i only for r = 0, otherwise a 
Dirac reduction is required. For r = by reduction of ( f4.29p we get the Benney system 



U \ n ( UU X ~\~ V x \ J tt - red 



t-2 



2 \ = n^dHx = 7r r n ed dH , (4.33) 



where 



. d x \ red _ ( 2d x d x u 
" _L " 1 ft. ' I ud r d^v + vd^ 



Hi= (u 2 v + v 2 ) dx, H = 

Jn Jq 



(4.34) 

uv dx. 



The next admissible reduction w = v = of (|4.28|) leads to 
2 - r 

u t2 _ r = y—^ uu x, (4.35) 

but for r = we lose the bi-Hamiltonian structure since there are no Dirac reductions 
with w = of ([4301 ). 

The case of k = 2. 

Example 4.8 Two field hierarchy: k = 2, r 6 Z \ {3}. 
The Lax operator is given by 

L = u 3 - r p 3 - r + vp 2 ~ r + p x ~\ (4.36) 

then we have 



4 — r 



u t4 _ r = 2 ^Zly ( 2 ( 3 " r )( 2r " 3)M 2 - r ^ + 2(3 - r)u 3 - r v x - (2 - r>V, 

+ (2-r) 2 (lnM) x t; 3 ), ( 4 - 37 ) 

4 - r 
2(3 -r) 2 

-2(1 -r) (3- r)u r ~ 2 vv x \ 



v u . r = ^ 2 ( r " !)(3 - r) 2 « s + (1 -r)(2 - r)(3 - r) M r "W 



For r = 2 we find 



u ) t2 \ u x + vv x 



u\ 2 u x v + uv x . 



which is again a Benney system with the known bi-Hamiltonian structure, this time 
reconstructed from formulae Q3.44 ) and ( [3.45| ). 
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Example 4.9 Two field hierarchy: k = 2, r G Z \ {2}. 

The Lax operator has the form ( |3.17| ) with N = 2 — r, m = r 

L = u 2 - r p 2 ~ r + vp l ~ r + p~ r , (4.39) 

Then we find 

^ = w^,L )r ~ ( : )^ = ( <v^r ) • < 4 - 4o) 

This Lax operator forms a proper submanifold as regards 7r_i only for r = 1, otherwise 
a Dirac reduction is required. Then for r = 1 



u \ uv 



v J. \ u, 

ti 



(4.41) 



it is again a Toda system with the known bi-Hamiltonian structure, this time recon- 
structed from formulae ( |3.44| ) and ( |3.46|) . 



Example 4.10 Three field hierarchy: k = 2, r G Z \ {2}. 

The Lax operator has the form ( |3.16| ) with N = 2 — r, m = r + l 

L = up 2 ' 1 ' + vp 1 ^' + wp~ r + p~ r ~\ (4.42) 

Then we find 

(« \ / (r - l)u x v + (2 - r)w x \ 
u = ru x w + (2 - r)uw x .(4.43) 

u;/^ V (1 + rK / 

This Lax operator forms a proper submanifold as regards n_i only for r = 1, in other 
cases a Dirac reduction is required. Then for r = 1 

uv x 

u x w + uw x ) = 7T_idffi = 7T T ed dH , (4.44) 
2u x 

where 







• 


'.-( 












ud x 





<9 x .u 














20* 



7T_ 



\ud x u ud x v \ u ^x w 
K ed = | u 3c w ^ + 9 x m + 2ud x 

\wd x u 2d x u + ud x —~wd x w + d x v + vd x 



(4.45) 



H x = I I (v 2 + 2ww) dx, H = I v dx. (4.46) 
2 in in 
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Example 4.11 Two field hierarchy: fc = 2,rGZ\{l}. 

The Lax operator has the form ( |3.17|) with N = 1 — r, m = r + 1 



L = u 1 ~ r p 1 ~ r + vp~ r + p r \ 



Then we find 

Lt 2 - T = 



2-r 
1,1-r 



L 



2-r 



ruu x v + ii D, 



1-r V (1 " r^u^u, 



(4.47) 



.(4.48) 



Let us consider the case of r = 0. To get 7r_i we have to make a Dirac reduction as the 
condition N > 2r — 3 > — m is violated. The simplest admissible Lax polynomial has 
the form 



L = up + v + wp 1 + zp 2 



and the Poisson operator reconstructed from (|3.44|) is 
/ 



7T-1 



2u 9 X - d x u \ 

ud x —v x 

(9 x u — d x w 

y 2d x u — u d x v x —w d x —d x z — zd x J 



Then, reduction of ( |4.5U| ) with constraints z — 0,w — 1 gives 



red 



7i_ 



u 2 d x 



d x u 2 

while the second Poisson operator, constructed from ( |3.45| ) takes the form 

u 2 d x u + u d x u 2 u 2 v x + u 2 d x v 
—u 2 v x + v d x u 2 2u d x u 

Fortunately, both Poisson operators are again differential. Hence 



7T, 



red 



*2 



u 2 v x 



^dH 1 = 7r r o ea dH 



red , 



where 



{u + v 2 ) dx, Hq 



v dx. 



Example 4.12 Three field hierarchy: k = 2, r 6 Z \ {1}. 

The Lax operator has the form ( |3.17| ) with iV = 1 — r, m = r + 2 



L = u 1 r p 1 r + vp ' + wp ' * + p 



-r-l 



-r-2 



Then we find 



2-r 
1,1-r 



L 



ruu x v + u v x 



1 — r 



ti- 



(1 — r)u l r {{\ + r)u x w + uw x ) 
(2 - r)(l - rju 1 ^ 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



(4.53) 



(4.54) 



(4.55) 



(4.56) 
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Let us consider the case for r = 0. Again condition N > 2r — 3 > — m is violated but 
reducing (|4.50[) with constraint z = 1 we get the first Poisson operator in the form 



red 



\ud x u — \u x d x 1 u x \uv x — \u x d x 1 v x —\uwd x + \u x w — \u x d x 1 w x 

* -\v x d~ x v x ud x + \v x w - \v x d~ l w x 

* * \w 2 d x + \d x w 2 - \w x d x 1 w x 



(4.57) 



where * denotes the elements that make the matrix skew-adjoint. The second Poisson 
operator calculated according to ( |3.45|) is 

« ed )n = \ u2 ( v ~ \ w2 ) 9 * + \l u ( v - \ w2 )* - u *( v - i w2 M lu * 

+ -d x u 2 {v - \v 2 ) + -u x d- l [u{v - ^w 2 ) x - u x (v - -w 2 )}, 
{K ed ) 12 = l u2wd x + \u{v - ^w 2 )v x + ^[u(v - ^w 2 ) x - u x (v - ^w 2 )]d^v x 

3 1 11 11 

(K ed ) 13 = (-jU 2 + gUw 3 )d x - -uu x + -w(2vu x - uv x - -w 2 u x + -uww x ) 

+ -[u(v - ^w 2 ) x - u x (v - ^w 2 )]d~ x w x - ^uvwd x 

I i 

+ ~^xd~ x \2u x - (vw) x + -(w 3 )J, 

( 7[ o ed ) 2 2 = \ ud * u + \[uw x -{v- ^w 2 )v x }d^v x + ^v x d x 1 [uw x - (v - ^w 2 )v x ], 

( 7r ed ) 23 = U ( V ~ \ w2 ) d * - \ \ U - W ( V - \ w2 )] V x - ^UWW X 

+ -Vxd^ftux - {vw) x + ^{w 3 ) x ] + ^[uw x - (v - ^W^Vxld^Wx, 

K^L = \[(vw 2 - 2uw)d x + d x (vw 2 - 2uw)} - ^(w% + d x w A ) 
+ -{2u - vw + ^w^) x d x l w x + ^w x d x 1 (2u - vw + ^w 3 ) x . 
Notice that both Poisson structures are nonlocal. Then, 







f u 2 v x \ 




1 -1 


uu x w + U 2 W X 









r r _ e fdHi = rt ed dH , 



(4.58) 



where 



Hi = / (uw 2 + v 2 w + 2uv) dx, H = (u + vw) dx. (4.59) 
Jo. Jn 
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